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Correlations of azimuthal angles observed at RHIC have gained great attention due to the prospect 
of identifying fluctuations of parity-odd regions in the field sector of QCD. Whereas the observable of 
interest related to parity fluctuations involves subtracting opposite-sign from same-sign correlations, 
the STAR collaboration reported the same-sign and opposite-sign correlations separately. It is shown 
here how momentum conservation combined with collective elliptic flow contributes significantly 
to this class of correlations, though not to the difference between the opposite- and same-sign 
observables. The effects are modeled with a crude simulation of a pion gas. Though the simulation 
reproduces the scale of the correlation, the centrality dependence is found to be sufficiently different 
in character to suggest additional considerations beyond those present in the pion gas simulation 
presented here. 

I. INTRODUCTION 

When modeling heavy-ion collisions, one usually ignores momentum conservation. Given the thousands of hadrons 
emitted in central collisions at RHIC (the Relativistic Heavy Ion Collider), such effects are of the order of a fraction 
of a percent. However, as will be illustrated here, high-quality data sets measured by the STAR collaboration 
have sufficient statistics and coverage to identify and analyze correlations at levels better than 10 -4 , enabling the 
identification and analysis of correlations due to momentum conservation. 

In the last year STAR has reported two-particle angular correlations of the form [TJ [2] , 

7 = (cos(0i + 2 )) = (cos 01 cos 02 - sin0! sin0 2 ), (1) 

where 0j refers to the azimuthal angle measured relative to the reaction plane. Correlations were reported for both 
same-sign, 7 SS , and opposite-sign, 7 os , pairs. The difference between the two correlations, 7 q bai = 7os — 7ss> was 
proposed to be a sign of parity fluctuations [3J 0]. The idea was that the parity-odd topological charge in the 
QCD sector would couple the large out-of-plane coherent magnetic fields from the passing ions to E a • B a in the 
QCD sector to create an out-of-plane electric field that would correlate out-of-plane charges of the same sign, or 
equivalently anti-correlate opposite-sign charges emitted out of plane. This would motivate a positive correlation for 
(sin 0i sin ss — (sin 0i sin 02) os - To reduce the contributions from resonances, the cos0icos0 2 correlations were 
then subtracted from the sin 0i sin 02 correlations, i.e., one would look for 7 q bai > 0. To better determine whether 
the observed angular correlations were evidence of large parity fluctuations, several more differential analyses have 
been performed: looking at the contributions to in-plane and out-of-plane separately [5,6 or looking a charge balance 
functions 0[S]. Such studies suggest that most of the correlation in 7 q b a i derives from another source. Upon further 
consideration it appears that most of the correlations in 7 q bai derives from charge conservation overlaid onto elliptic 
flow [7]. In fact, the more differential correlations encoded in charge balance functions [5] were even reproduced at 
the quantitative level by overlaying local charge conservation onto a thermal blast-wave model which was tuned to 
single-particle spectra and elliptic flow. 

The effects of charge conservation overlaid onto elliptic flow only addressed 7 q bai and are irrelevant for explaining 
why the same-sign correlation, 7 SS , was fairly large and negative. Like the values for 7 q bai> the magnitude of the 
correlation was found to be of the order V2/M, where M is the multiplicity and V2 — (cos 20) is the elliptic flow. It 
was argued in [10] and that such correlations might ensue from momentum conservation combined with elliptic 
flow. If for each particle emitted at angle 0i, there was a particle balancing the transverse momentum at angle 
02 = 0i — 7r, one would find correlations, 7 ~ —V2/M. 

The goals of this paper are: (1) to expand the discussions of momentum conservation in [101 and to discern 
what combination of 7 SS and 7 os , and with what weighting in momentum, would best isolate the effects of momentum 
conservation. (2) to demonstrate the effects of momentum conservation with a simulation, compare to experiment, and 
see whether the observed correlations are consistent with expectations from momentum conservation. (3) Investigate 
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whether correlations of the class described above can be used to discern fundamental properties of the matter created 
at RHIC, e.g., the viscosity. The next sections presents the principles involved in observing momentum-conservation 
driven correlations with a focus on the momentum sum rule which drives the phenomenology. Section |TTl| discuses the 
role of viscosity in spreading the correlation in relative rapidity. The pion cascade model is presented in Sec. |IV| and 
compared to experimental results from STAR, while a summary and outlook are given in the final section. 



II. THE MOMENTUM SUM RULE 



In the absence of competing correlations, such as charge conservation, the correlation (cos(0i + ^2)), averaged 
over all pairs regardless of charge, would be well suited to study momentum correlations. Even better would be 
the correlation, (p x iPx2 ~ PyiPy2) / (p 2 ) , which is similar but is determined more directly by momentum conservation. 
Unfortunately, other correlations come into play which makes the choice of observable tricky. For instance, one might 
consider same-sign correlations as a way to avoid the effects of charge conservation. However, when one sees a positive 
particle at 4>i, there is an enhancement for having a negative particle at a similar angle and with similar momentum. 
The momentum of the negative particle also needs to be balanced, and the balancing momentum could be shared 
by other positive particles. One must also consider the fraction of the balancing momentum carried by unobserved 
particles. For instance, in the STAR analysis neutral particles are not measured and even with the large acceptance 
and relatively high efficiency a good fraction of the balancing momentum is outside detection. 

To better clarify the issues above, the theory discussion is split into parts. In the first subsection the situation 
will be considered for a system with a single species. The momentum sum rule will be derived for perfect acceptance 
and efficiency, with a discussion of how it would be modified for less than perfect detectors. The more realistic case, 
accounting for undetected neutral particles and the effects of charge conservation, is then considered in the following 
subsection. The relations found here are applied and tested with a simple model in Sec. |IV| 



A. The Idealized Case with No Other Correlations 



For illustrative purposes we consider the case of M particles emitted under the constraint of momentum conservation. 
Rather than calculating 7 as defined in Eq. ([lj, we consider the p t weighted correlation, 

, = (PxlPx2 ~ PylPy2) (t) , 
7 " <P?) ■ ^ 

If all the particles had the same pt, this would be no different than 7, which is defined solely through the angles. 
Although the angular correlations are stronger for higher p tl 7' surprisingly turned out to be several tens of percent 
smaller than 7 in the simulations of the next section. This is due to the denominator (p 2 ), which weights high- 
momentum particles even more strongly in the denominator than the extra pt weightings in the numerator. If the 
denominator had been defined as (pt) 2 rather than (pf), 7' would have been larger. The particular form in the 
definition of Eq. ([2]) is chosen to simplify the momentum sum rule shown below. Since this form is most directly 
related to the sum rule, it exhibits less model dependence than other choices. 

To derive the momentum sum rule, one expands the definition of the averages assuming for the moment that all 
particles are measured, 

, 1 Ysij&j \PxiPxj — PyiPyj) , . 

7 = SI (3) 

1 Si j (PxiPxj ~ PyiPyj) 



M-i EiPi 
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1 Sj=j {PxiPxj PyiPyj) 



M-l I>? 

2 2 
^xi Pyi 



1 s l,i>i, 1 J2iPti cos2 <l>i 



M-l M-l E,p? ' 

Momentum conservation, ~Y^iPxi — 0, was applied in the third line. The final term resembles an expression for i! 2 
weighted by p\ . Assuming M is large one can derive the momentum sum rule, 



(4) 
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The sum-rule above is based on a perfect acceptance. Even for the STAR TPC detector, which has a broad 
acceptance of approximately ±1 units of rapidity, only about one third of the emission is addressed. Since the 
averages (• • • ) are correlations, the sum rule is not affected by a less than perfect efficiency if the efficiency is uniform. 
However, if the experiment has a limited acceptance, the sum rule can be significantly compromised. First of all, 
one must state whether M refers to the multiplicity inside the acceptance, to the total multiplicity, or something 
in between. For instance, the STAR experiment has good coverage for pseudo-rapidities between -1 and 1 and for 
charged particles with transverse momenta above ^200 MeV/c. One might choose M as the multiplicity of particles 
with p t > 200 MeV/c, or one might use the usually reported measure of the multiplicity which assumes no pt cuts. 
If M refers to the multiplicity within the acceptance, the sum-rule is modified by a factor f p , where f p is defined as 
the fraction of balanced momentum which falls within the acceptance, 



Here, fp < 1, and the sum rule becomes, 



Part of the reduction, fp, comes from balancing momentum that is carried by particles with rapidities outside the 
acceptance. This can be estimated by looking at 7' as a function of relative rapidity and then convoluting that 
with the rapidity distribution. Given that the efficiency inside the acceptance will vary, the surest way to proceed 
is to model the correlations and simulate their measurement. Unfortunately, either strategy introduces some model 
dependence, which is unavoidable once acceptances become complicated. 

Although transverse momentum is conserved globally, the balancing momentum for a specific particle might be, 
or might not be, spread over the entire collision volume. Whereas STAR's analysis cover approximately ±1 units 
of rapidity, emission is spread over a half dozen units of rapidity. In the thermal calculations by [llj . transverse 
momentum conservation was enforced globally. If one were to enforce conservation locally, the question would be 
"How local?". Clearly, addressing this question is central to estimating the fraction fp in Eq. ([5]). During the initial 
interactions of the incoming nuclei, hard scatterings and fragmentation processes spread transverse momentum over 
a wide range of rapidity, whereas according to the commonly accepted picture of a RHIC collision, the subsequent 
processes are more local in nature, as the momentum transfers from subsequent collisions are shared amongst particles 
with similar rapidities. Further, it is these collisions that build up elliptic flow. Here, we will show how the lack of 
local momentum balance in the initial scattering does not affect the locality of Eq. ^ as long as the particles created 
in the initial state have their momenta distributed isotropically in azimuthal angle. 

First, one considers a distribution of particles with initial momentum ki, whose final momentum is pi and whose 
net momentum transfer is ft, 




(5) 



Pi = h + q % . 



(7) 



We will assume that the vectors fc, are isotropic, i.e., 




(8) 



and that the momentum transfers conserve momentum within the neighborhood 




(9) 



Here, fp is again the fraction of the balancing momentum exchanges, ft, that will be found within the acceptance. In 
the limit that the acceptance covers much more phase space than what is subtended by the diffusion of momentum 
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from the transfers g,, fp — ¥ 1, even if the collision then covers many more units. Calculating the azimuthal correlation, 
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M(M-l) 



*2kyiQyj -\- QxiQxj QyiQyj) 

M(M — 1) (^kxiQxj — 2k yi q y j + QxiQxj — QyiQyj) 



-fp 



M(M - 1) 

-fp 



M 



(2k x q x -2k v q y + q x -q y ), 



where the isotropy of k and local conservation of q were both exploited. For the purposes of this derivation we consider 
a class of events with fixed M, but will take the limit of large M . 
One can also write an expression for the numerator for v' 2 , 

«f Um) = (k 2 x-kl + 2k x qx-2k y q y + ql-q 2 y ) (11) 
= (2k x q x - 2k v q y + q% - q%), 

where the isotropy of k was again used. 

The equivalence of the the last expressions in Eq.s (10) and (11) confirms Eq. ^5J. If the momentum transfers 



responsible for elliptic flow build up during the first few fm/c of the reaction (but not in the initial state), one would 
expect the extent of 7' as a function of relativity to depend on the diffusion of transverse momentum along the 
longitudinal direction, and thus perhaps on the shear viscosity. This will be explored in the next section. 

B. Including the Effects of Charge Conservation 

Charge conservation affects both same-sign and opposite-sign correlations. It is straight-forward to understand 
how opposite-sign correlations become manifest. Here, we first describe how charge conservation also affects same- 
sign correlations. Due to local charge conservation, each charge (aside from those carried by the incoming nuclei) is 
accompanied by an extra charge of the opposite sign. Since charge conservation is local, the two balancing charges 
will remain close to one another in coordinate space. The two balancing charges then flow together and are positively 
correlated in <j). If one triggers on a specific positive particle, one expects momentum conservation to provide negative 
correlations in azimuthal angle, and that the correlation would be spread out roughly equally among positive, negative 
and neutral particles. The same-sign correlations would be expected to account for one third of the momentum balance. 
However, that balance is increased by the fact that the original positive particle is accompanied by a negative particle 
with similar momentum. This accompanying particle's momentum must also be balanced, and if it has the same 
momentum as the original triggered particle, the momentum balance one would find amongst other positive particles 
would double, i.e., one would expect nearly 2/3 of the original momentum of a triggered positive particle to be 
balanced by other positive particles. 

The first goal of this section is to find a combination of 7g S and "f' os that leads to a momentum sum-rule despite the 
correlation inherent due to each positive being balanced by a negative charge, and the complication due to roughly 
one third of the particles being neutral, and in the case of STAR, undetected. To proceed we consider a system with 
multiplicities M± for both positive and negative particles and multiplicity Mq of neutrals. Assuming that charged 
particles correlate equally with neutrals as they do with particles of the same charge, one can write the various 
correlations as: 
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'l2ie+,j<e+,i^j^P xi P x 3 PyiPyj) + ^2ie-.je-,i^j^P xi P x J PyiPyj) 
2M±(M ± -l)(p?) 

+ {PxiPxj — PyiPyj) 



7ss — r, M I n ,r TT7Z2\ ' 



Ml(p? 



i _ 12ie+,jeo(.P x iP x j PyiPyj) _ , _ , 
7+0 ~ M ± M {p 2 t ) ~ 7 -° ~ 7ss 
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Using the symmetries listed above and assuming the multiplicities are large numbers, 

xi -P 2 yi) 



(M ± +MoH s + M ±7 > s = - (13) 



2M± (p: 



(M± + M )7b S + M±j' os = -v' 2 . 

If M refers to the net multiplicity of charged particles, M = 2M± , and if /q is the fraction of particles that are neutral, 
Mo = fo(M + Mo), one express the result as: 

1A+/o \;+k s = -i d4) 



2 V 1 -/o/ 2 ,os M 
For one third of the particles being neutral, /o = 1/3, 

7ss + ^ = ~M- (15) 

Since this linear combination of same-sign and opposite-sign contributions seems most connected to the momentum 
sum rule, we define the quantities 

/ / i 7os / 1 n \ 

7 P bai = 7 SS + -y, (16) 

. 7os 

7pbal = 7ss + -y ■ 

By explaining both 7 q bab which focuses on correlations from charge conservation, and 7 p bai which is most related to 
momentum conservation, one would reproduce both the same- and opposite-sign correlations. The advantage of this 
choice is that the individual quantities can best be interpreted in terms of a single effect. 

Following the lines from the previous subsection, one can also account for lack of perfect acceptance, at which point 
one finds the following final expressions, 

HW^J 788 2 7 ° S = ~ fP M' ( ] 

+ — = -fp V f r for /o = 1/3. 

Here, the multiplicity M is the actual measured multiplicity, and fp is the fraction of balancing momentum observed 
given the limited acceptance and efficiency For a constant efficiency in a limited acceptance, efficiency reduces both 
fp and M by the same factor, as expected since correlations should be unchanged by a constant efficiency. Thus, this 
expression also works if M is the total multiplicity inside the acceptance assuming perfect efficiency, as long as fp 
refers to the fraction of momentum inside the acceptance, not just what is measured. For complicated efficiencies and 
acceptances, the model independence becomes somewhat compromised. However, even in that case the expression 
suggests that the combination j' as + r y' os /2 is a good choice for an observable that would best isolate the effects of 
momentum conservation. Additionally, one must consider the qualifier that analyses of i>2 yields a variety of results, 
differing by 10-20%, depending on how the analysis was performed. For the purposes of the sum rule, the value for 
V2 should include some "non-flow" contributions, such as those arising from resonances. 



This is a good point for summarizing the assumptions that led to the derivation of Eq. (17 1. First, it was assumed 
that the momentum transfers responsible for elliptic flow are locally balanced, i.e., not spread over many units of 
rapidity. Secondly, it was assumed that an unmeasured neutral particle would balance the momentum of a given 



charged particle identically as another charged particle of the same sign. If observations are not in line with Eq. (17 1, 
one should question these two assumptions. 

III. MOMENTUM DIFFUSION AND VISCOSITY 

It was shown in |12H14j that transverse momentum diffuses longitudinally with the viscosity playing a role in 
determining the diffusion constant. Restating the arguments, one can consider a boost-invariant system with a 
momentum per rapidity slice, Px{v) or Py^fi^ wnere 



Px{v, T )^ T dxdyT 0x (x,y,rj,T), (18) 



G 



where T a p is the stress-energy tensor as viewed by an observer moving with rapidity, y — r/. Here, 77 and r play role 
of the longitudinal position and the proper time in Bjorken coordinates, 

z = rsinh?7, £ = tcos1i?7. (19) 

Momentum conservation, O t Tq x + diT vx = 0, yields 

9 T P X (7?, r) = -ify J dxdy f xz . (20) 

Assuming the Navier-Stokes equation, 

f xz = - r ^-d rj v x , v x = e ^°p ) (21) 
where ry s is the shear viscosity, one then finds 

d r P x (v,r) = Tj^j *2 9 V P *(V, r). (22) 
For small z, (l/r)d n — d z , so the diffusion coefficient becomes 

D =j?Te- (23) 

For the correlations in a boost-invariant system, we define the correlation function for the interval < r\ < H , with 
the rapidty range H — > 00, 

1 



IUA77) = - J dr h drk((P x {r} 1 ,T)P:{r}2,T))).8(ri 1 -rh-&r]) (24) 



where the averaging, ((• • • )), denotes an average over events. If one subtracts the i = j terms, one could identify the 
sum as the numerator used for calculating 7p ba j- This motivates the definition of 11^,, 

II^Ar?) = U X (A V ) i /<V P M W? ? ) (25) 



H 



Using Eq. (22 1, the evolution of W! x is 



2D 

d T IL' x (Ar l ) = ^d 2 Il'(Ar 1 ), Ai^O. (26) 



The condition that one must stay away from Arj = comes from the diffusion equation not being valid for scales 
so small that the individual particles correlations with themselves does not appear. To include Arj — 0, one can 
consider the effects of collisions on II' defined in Eq. (25). Comparing the instants immediately before and after a 
collision between particles i and j, the contribution to II' involving other particles than those involved in the specific 
collision, k 7^ i, k ^ j , does not change because the total momentum carried by i and j does not change. However, 
the contribution to n' involving the colliding pair, (z, j), does change suddenly. To calculate the change, consider the 
fact that the net momentum does not change, i.e., 

A( Plx + Pjx ) 2 = 0, (27) 

One can then identify the change in n' as 

AU'^Arj) = ~A(2p ix p jx )5(A v ) (28) 



~A{p 2 ix +p 2 jx )6(A v ), 
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where it is assumed that the two scattering particles have the same coordinate r\. Similar expressions can be derived 
for any m —> n process. The rate of change of II' due to the instantaneous collisions thus has two pieces. First, 
the motion of the particles in between collisions is responsible for the diffusive term. A second term involving the 
instantaneous changes due to collisions of two particles at the same rapidity can be identified as the collision rate 
multiplied by the average change of p% per collision, 

d T IL' x (Ar,) = ^d 2 v K(A v ) - Ai ( ) ) *(Atj). (29) 



The second term in Eq. ( 29 1 is a source term for the diffusion. The quantity IL' X — II' is proportional to the 
numerator used to define 7', as m Eq. ([2]), if it were binned for pairs of a specific relative momentum. Thus, 

7 (A,, r) = pi. - II,) = , (30) 

where {pf) is defined for large asymptotic times. 

Since 7' is proportional to II^, — Tl' y , which obey diffusion equations, so does 7', 

d Tl '(A V ,r) - ^7'(A*7,t) = -^(A^v'^r). (31) 

The source term for the diffusion represents the rate at which v' 2 rises, defined in Eq. Q , with all the strength located 
at Aij — 0. The source term thus maintains the momentum sum rule for 7' derived in Eq. Q. 

To gain an initial understanding of the scale one would expect the diffusion to extend, one can consider a boost- 
invariant system of multiplity, dM/dr], in area A, of particles colliding with fixed cross sesction a. The diffusion 
constant increases with time due to the fall in density. For massless particles, one can crudely express the viscosity 
in terms of the cross sections and density as, 

(32) 

5 na 

If the density is (dM / drf) / {At) and if one considers times for which the source term is zero, the diffusion equation 
becomes 



A solution to the equation is 



k exp-(A,)V4^ 

- 1 A- 2A < h 

s — m r, a 



{dM/drj)a' 

where the transverse area is A. The width of a diffusive cloud that beg an as a delta function at T c }}jp^j C is, 



A/5 

ln(r/T ni p tic). (35) 



' y (dM/drj)a 

Thus, Tdiiptic would represent a typical time for which elliptic flow is generated. This tends to be early in the collision, 
perhaps a few fm/c after thermalization. For central collisions, the transverse area is a few tens of square fm, and the 
multiplicity can be several hundreds. If the characteristic cross sections are a few tens of mb, and if the ratio r/r e iiiptic 
is of order ten, the spread should be roughly a half unit of rapidity. This is the same scale as the spread due to the 
final thermal motion. The breakup conditions and flow can be well understood by either blast-wave analyses, or by 
fitting to hydrodynamic or to hybrid hydrodynamic/cascade models. This would allow isolation and determination of 
as a function of centrality. If taken in concert with a detailed model of the evolution that describes when elliptic 
flow was generated, this might provide insight into the viscosity. 

It should be emphasized that the simple calculations of this subsection, based on a boost-invariant model with no 
transverse flow, and with an assumption of a fixed number of massless particles with fixed cross sections, certainly 
cannot be seriously applied to extract properties such as the viscosity from data. However, this simple picture provides 
insight into the role of viscosity and the importance of understanding the origins of elliptic flow. 



IV. PION CASCADE 



In order to investigate the ideas put forth in the previous sections, and in order to see whether the ideas might 
ultimately explain measurements at RHIC, a simple model was constructed in which a gas of pions evolved according 
to simple s-wave cross scatterings with a fixed cross section. The initial momenta were assigned according to a thermal 
distribution with a uniform temperature. The space time coordinates were assigned according to the distribution, 



with the collective velocity in the z-direction given by v z = z/t. These are known as Bjorken coordinates, which 
are invariant to boosts in the longitudinal direction, similar to a flat-space Hubble expansion, in the limit r\Q — > oo. 
After picking source points according to the distribution above, two-particles were generated according to the thermal 
distribution for each point, then boosted according to the longitudinal collective velocity. Two thirds of the pairs 
were chosen to be ir + ir~ pairs, and one third as 7r°7r° pairs. This enforced local charge conservation at the extreme 
in the initial state, which then relaxed as the system evolved. The collective nature, short mean free paths, led to 
strong correlations between the balancing particles in both azimuthal angle (j> and in rapidity. The evolution involved 
straight-line trajectories punctuated by s-wave collisions with fixed cross sections that were independent of isospin. All 
the simulations used values of Tq = 250 MeV for the initial temperature, tq = 1.0 fm/c for the initial time and t\q = 2. 
The fixed cross section a was set at 15 or 30 mb to test sensitivity of the final-state correlations to diffusion. Runs 
were made with 1000 particles and with initial transverse radii of both R x = 1.5, R y = 3.0 fm, and R x = 1.8, R y = 2.5 
fm. These two geometries give the same central densities, thus making it possible to test the scaling with anisotropies. 
Additionally, runs were made with a smaller system, 250 particles, with radii of R x — 0.9, R y = 1.25 fm, in order 
to vary the size while keeping the density fixed. Due to the high initial pion density, ~ 1 pions/fm 3 , and the large 
cross sections, a = 30 mb = 3 fm 2 , momentum transfer can diffuse at super-luminous speeds [151 116) . To prevent 
such unphysical behavior the particles were oversampled by a factor iV samp i c combined with a reduction of the cross 
section by a factor 1/N sa . mp i e . In the limit -/V samp i c — > oo, super-luminous behavior would disappear while maintaing 
roughly the same mean free path as for jV samp i c = 1. For the 15 mb case, it was found that -/V samp i e = 4 was sufficient 
to approach the N^^i^^ limit, while iV samp i c = 8 was used for the a — 30 mb calculations. Analyses were made 
for two acceptances: a perfect acceptance and one that crudely mimics STAR's acceptance. The STAR acceptance 
was modeled by considering only those particles with pseudo-rapidities between -1 and 1, and that had transverse 
momenta between 200 MeV/c and 2 GeV/c. 

Results for M^ pbsd /v' 2 and A/7 p b a i/«2 ar e shown in Table |lj for both the case of perfect acceptance and for STAR's 
acceptance, for both anisotropies, for both cross sections, and for the scaled-down size. The multiplicity M is 
the average number of charged particles with pseudo-rapidities between -1 and 1. The statistical accuracy of the 
quantities are of the order of 5-10%. Calculations with smaller i>2, or with smaller acceptances, tended to have larger 
uncertainties. For that reason, runs were not made with smaller anisotropies. 

The momentum sum rule applies only to the case of full acceptance, and only for the primed quantities. The 
results in Table |l] show that finite acceptance reduces the ratios, Ip^M/v^, by approximately 20%. This reduction 
implies that for each particle within the specified rapidity range, roughly 80% of the balancing transverse momentum 
exchanges responsible for elliptic flow can also be found within the rapidity range. The reduction appears fairly 
uniform for all the cases studied. Unfortunately, STAR did not report the primed (p t weighted) quantities in [2], but 
instead reported unprimed quantities. Fortunately, the unprimed quantities in Table [I] are also fairly constant. The 
steadiness of the values would lead one to expect values for 7 p bai-^/^2 to be near 1.3, and to vary rather little with 
centrality. Since the ratio depends on the fraction of balancing momentum found within the acceptance, and since 
that fraction might well depend on centrality, it would not be surprising if the ratio varied by tens of percent. 



Experimental results for the ratio, r y p b&\M/v2, are displayed in Fig. [Tj The expectations described above are not 

well validated by the experimental observations. Although the sign ancToverall magnitude are of the right order, the 
experimental ratio increased by a factor of three from peripheral to central collisions, whereas the ratios extracted 
from model varied rather little from 1.3, roughly independent of changes in the initial anisotropy, the cross-section, 
or the system size. This discrepancy between the model and experiment at the qualitative level suggest a significant 
shortcoming in the model. 

Aside from the large centrality dependence, a second surprise is in the size of the ratio for the most central collisions. 
For ratios much smaller than 1.3, a simple explanation would involve a larger fraction of the balancing momentum 
being carried outside the acceptance. For ratios much higher than 1.3, one would look for reasons to expect a larger 
fraction of the momentum being balanced by particles inside the acceptance than what was occurring in the pion 
cascade. However, for the pion cascade ss 80% of the momentum sum rule was being realized, which gives little room 




(36) 



T = 



\lt 2 - z 1 , rj = tanh _1 (z/i), 
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acceptance 


a (mb) 


R x , R y (im) 


«2 


7pbalM/"2 


7pbalW«2 


FULL 


15 


1.8, 2.5 


0.0677 


-1.901 


-0.985 


FULL 


30 


1.8, 2.5 


0.0776 


-2.034 


-1.060 


FULL 


15 


1.5, 3 


0.137 


-2.075 


-1.007 


FULL 


30 


1.5, 3 


0.158 


-2.051 


-1.002 


FULL 


15 


0.9, 1.25 


0.0503 


-1.798 


-0.982 


FULL 


30 


0.9, 1.25 


0.0562 


-1.959 


-0.983 


STAR 


15 


1.8, 2.5 


0.0910 


-1.324 


-0.817 


STAR 


30 


1.8, 2.5 


0.0997 


-1.307 


-0.779 


STAR 


15 


1.5, 3 


0.186 


-1.282 


-0.789 


STAR 


30 


1.5, 3 


0.204 


-1.407 


-0.847 


STAR 


15 


0.9, 1.25 


0.0667 


-1.1354 


-0.763 


STAR 


30 


0.9, 1.25 


0.0697 


-1.2854 


-0.833 



TABLE I: The results for i>2 and the scaled correlations Mf p ha,i/v2 and M^' pbal /v2 are shown for different acceptances, cross 
sections, initial radii and sizes. The scaling (multi plied by the multiplicity /t^) was chosen so that 7pb a i-^/' !; 2 would be unity 
for full acceptance according to the sum rule, Eq. (15 1. For the calculations with STAR acceptance, the values of 7' hovered 
around 0.8, indicating that for any given observed particle about 80% of the balancing momentum responsible for elliptic flow 
would be found within — 1 < 77 < 1. The scaled values of 7 were consistently nearly twice as large as the values for 7'. Given 
the statistical errors, of the order of 10% using 40,000 events, it is difficult to see any sensitivity of the scaled moments to 
doubling the cross section, doubling the system size, or roughly doubling the anisotropy. Given the consistency of the 7pb a i an d 
7 P bai values, these calculations would lead one to expect similar behavior in STAR analysis, i.e., the scaled values of 7 would 
stay roughly constant with centrality, and within a few tens of percent of 1.3. 
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FIG. 1: (color online) Results from STAR are shown for the ratios M7 P b a i/w2 and for M7 q b a i/«2. The ratios M7 q b a i/«2, 
which are related to momentum conservation, vary by a factor of 3 from the most peripheral to the most central collisions. In 
contrast, the values found with the cascade model (horizontal lines) vary only by 10%. The parameters for each of these values 
is described in detail in Table [T] with the variation being statistically consistent with being constant, even though geometries, 
sizes and cross sections were all varied by a factor of two for different runs. The ratio M7 q b a i/i>2 is identified with charge 
conservation. Despite the fact that the model poorly matches spectra or flow variables, the values found with the model were 
within a few tens of percent of the ratio found experimentally. 
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0.5 1 1.5 2 

AT| 

FIG. 2: (color online) Values of the scaled correlation, Mf p ha,i/v2, as determined by STAR, are shown as a function of relative 
pseudo-rapidity. Results are displayed for two centralities, 10-30% (red filled stars), and 30-50% (green open stars). The 
experimental correlations appear broader for the more peripheral collisions. Also shown are results from four of the cascade 
simulations: (blue squares) default results with initial sizes of R x = 1.8, R y = 2.5 fm and 15 mb cross sections, (purple circles) 
same as default but with double the cross section, (cyan diamonds) same as default but with half the size for R x and R y and 
with 1/4 the multiplicity, (brown downward triangles) same as default but with a larger anisotropy, R x — 1.5, R y = 3.0 fm. The 
model roughly reproduces the size and extent of the correlations. However, the cascade calculations show little sensitivity to 
any of the changes, at least within the statistical accuracy of the calculations. Although the statistical significance is marginal, 
it appears that the fall of 7 P bai for peripheral collisions seen in STAR's results, plotted in Fig. [l] derives from an increased 
broadening. 

The scaled moments M7 q b a i/»2, which are associated with charge conservaton, are also displayed. They also have a similar 
width and height as the experiment. Unlike the scaled values of 7 P b a i, the width for these correlations are sensitive to the cross 
section in the cascade model. 

for overshooting the ratios from the pion cascade. Another possibility would be that a larger fraction of the balancing 
momentum was being carried by neutrals in the real collisions as compared to the presumed 1/3 of the pion cascade. 
Yet another possibility might involve the ratio between 7 and 7' in Table This ratio might change in realistic 
models, though it is difficult to foresee by what amount. The better way to handle this would for the experiment to 
evaluate 7 pba i instead, as it would reduce the model dependence involved in interpreting the result. The experimental 
ratio in Fig. [l]was approximately 40% higher for the most central collisions than in any of the cascade calculations. 
Systematic errors for both i>2 and for the angular correlations were well over 10%, so the overshoot might disappear 
with a more detailed analysis or with higher statistics, but for the moment seems puzzling. 
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STAR has also presented more differential results for the moments displayed in Fig. [T] In Fig. [2j results are 
shown for the scaled moments, 7 P bai-^A?2 binned as a function of pseudo-relative rapidity. Analyses from STAR 
for centralities of 10-30% and for 30-50% are plotted along side pion cascade results for the initial anisotropy R x = 
1.8, R y = 2.5 fm, and for two cross sections. As can be seen in Fig. [2j the width of the correlation in reasonably well 
reproduced by the cascade. This centrality range corresponds to the values for which the net moments were reasonably 
well reproduced by the cascade as shown in Fig. [T] For the experimental results, it does appear that the correlations 
are more spread out in relative rapidity for less central collisions, but it would be useful to observe the correlations 
for both the most central collisions, and for the more peripheral ones, where the net moments are not in line with the 
expectations of the pion cascade. Performing the more differential analysis in relative pseudo-rapidity for a broader 
range of centralities could validate whether the dependence of 7 P bai with centrality indeed comes from a broadening 
of the correlation for peripheral collisions. Since the sum-rule seems rather robust from a theoretical perspective, it 
seems likely that the weakening comes from a broadening, but the experimental evidence for this statement needs to 
be strengthened. 

For the moment, we will assume the preliminary interpretation above is true: that the momentum balance re- 
sponsible for the generation of v 2 is local for central collisions, but much broader for peripheral collisions. This 
behavior is not matched in the pion cascade model, and given that the effects of momentum diffusion increase only 
logarithmically with time as shown in Sec. [Til it would seem most unlikely for momentum diffusion to substantially 
increase the width of the correlations in Fig. 2] for peripheral collisions, unless the processes responsible for elliptic 
flow occurred within the first few tenths of a fm/c. At such early times it would be difficult for elliptic flow to be 
generated hydrodynamically, even if hydrodynamics were valid, since the characteristic times for generating flow is set 
by the transverse density gradients. Another possibility is that the assumption that the source for 7' in the diffusion 
equation, Eq. |31[ is not really a delta function. In the summary we discuss the possibility that transverse momentum 
conservation related to elliptic flow is not purely local in rapidity. 



V. CHARGE BALANCE CORRELATIONS FROM THE PION CASCADE 



The correlation, 7 q bai = 7ss — 7os> is largely dominated by charge balance effects. Such correlations are driven by 
local charge conservation overlaid onto collective flow. These correlations can be well reproduced with a blast-wave 
calculation tuned to reproduce elliptic flow and spectra if the thermally generated particles are always produced in 
small neutral ensembles instead of singly [7J. Not only do such calculations give reasonable values of 7 q b a i, but they 
are remarkably successful at matching more differential observables such as charge balance functions. Since the pion 
cascade model presented here is not well tuned to reproduce flow observables or spectra, one should not take the 
success, or lack thereof, in reproducing data too seriously. Nonetheless, for the sake of completeness, we also show 
experimental results for 7 q bai in Fig. [T] along with the values coming from the pion cascade model, which are also listed 
in Table In] Since non-zero 7 q bai also requires elliptic flow, and since the effects fall inversely with the multiplicity, 
the scalcdvalues, r y q baiM/v2, are presented as opposed to 7 q bai itself. The experimental values of 7 q bai-M jv-i in Fig. W 
vary between 1.0 and 1.4, whereas the result for the pion cascade tend to be modestly lower. Unlike the scaled values 
of Tpbai, there are no sum-rules that apply for 7 q bai- As expected, the scaled values of 7 q bai are higher for cascade 
calculations with higher cross sections, or for larger systems. In such cases, balancing charges will by more strongly 
focused into the same direction. The cascade calculation appears to modestly under-predict the data, perhaps by 
~ 20%. This might suggest that in the real collision that many balancing charges are made later in the reaction than 1 
fm/c, which would lead to the charges being even more tightly correlated. However, that conclusion is premature given 
the numerous shortcomings of the model. The only firm conclusion to be made here is that the sign and magnitude 
of 7 q bai from the model and from data are fairly similar. 

Since the cascade model employs isospin-independent cross section, the charge balance correlations can be calculated 
by comparing only a particle with its balancing partner, thus avoiding the noisy statistical subtraction involving charges 
from different balancing pairs. Thus, the statistical accuracy for the calculation of 7 q bai is much better than that for 
7pbai- More quantitative insight into the data can better be ascertained with models that match both the spectra, 
yields and flow variables, such as what was studied in [7J. 



VI. DISCUSSION, SUMMARY AND OUTLOOK 

As predicated in the previous section, if observations are not in line with the predictions of the model, one should 
reconsider the two basic assumptions used to derive the sum rule, 7 p b a i-^/ w 2 ~~ * 1j f° r a broad acceptance, or the 
expectation that 7 P bai ~ 1.3 for STAR's acceptance and only weakly dependent on centrality: 
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acceptance 


a (mb) 


R x , R y (fm) 


«2 


7qbalM/u 2 


FULL 


15 


1.8, 2.5 


0.0677 


1.351 


FULL 


30 


1.8, 2.5 


0.0776 


1.546 


FULL 


15 


1.5, 3 


0.137 


1.326 


FULL 


30 


1.5, 3 


0.158 


1.524 


FULL 


15 


0.9, 1.25 


0.0503 


0.9816 


FULL 


30 


0.9, 1.25 


0.0562 


1.220 


STAR 


15 


1.8, 2.5 


0.0910 


1.056 


STAR 


30 


1.8, 2.5 


0.0997 


1.257 


STAR 


15 


1.5, 3 


0.186 


1.035 


STAR 


30 


1.5, 3 


0.204 


1.235 


STAR 


15 


0.9, 1.25 


0.0667 


0.844 


STAR 


30 


0.9, 1.25 


0.0697 


1.051 



TABLE II: Results from the pion cascade model for the scaled correlations 7 q bai are shown for different acceptances, cross 
sections and initial radii. These correlations are driven by charge conservation. 



1. It was assumed that for a charged particle, finding the balancing momentum related to elliptic flow in other 
particles of the same charge would be identical to finding the balancing momentum in neutral particles. This 
assumption was explicitly enforced in the pion cascade, but could be different in a more realistic picture. For 
instance, for positive pions coming from p + decays, one would have a contribution in the same direction for 
neutral pions, which would not exist between the positive pion and other positive pions. The negative correlation 
with same-sign pions would then be stronger to compensate. 

2. The momentum transfers that build elliptic flow are locally balanced, i.e., if a particle scatters into the reaction 
plane, the balancing momentum is likely to be found nearby in rapidity. For the momentum to be balanced 
by particles far removed in rapidity, the transfers need to have occurred early in time, or perhaps even during 
the passing of the initial nuclei. The speed for which fluctuations of the transverse momentum can traverse the 
matter in the longitudinal direction is limited by the viscosity in a hydrodynamic context [T2TI14) . 

Violations of (1) would lead to stronger correlations than those predicted in the model. However, it is difficult 
to expect that by correcting the model with more realistic resonance contributions one might find corrections more 
than a few tens of percent. Thus, such considerations might help explain why STAR's observations overshot the pion 
gas by «40% for the most central collisions in Fig. [TJ but are unlikely candidates for explaining the factor of three 
variation in Mj q b a \/v2 observed by STAR and displayed in Fig. [I] 

Violations of (2) should lead to smaller values for 7 P b a i for a unite acceptance. For the pion cascade calculations, 
the particles were initialized with zero transverse flow and with symmetric stress energy tensors. The width of the 
resulting correlations was remarkably insensitive to the cross section employed in the cascade. A shorter mean free 
path gives a smaller viscosity and thus a smaller diffusion of transverse momentum fluctuations, but also keeps the 
matter together longer. The surest way to increase the width of the correlations in relative pseudo-rapidity would be 
to begin the cascade earlier. From the discussion of the diffusion of momentum in Sec. |III[ an earlier start time would 
broaden the final width in the coordinate r\ logarithmically with the starting time. The solutions to the diffusion 
equation would sugggest that to double the width, one would need to start the cascade at times near 0.1 fm/c. These 
would be well below the times that would be reasonable for modeling the evolution as a cascade, even with partons. 
Further, one would need to explain how elliptic flow was being generated at such short times as the characteristic 
time should be set by the transverse size of the collision region divided by the speed of sound. 

It is difficult to judge how much the ratio 7pbai/7 p b a i ml ght change with a more realistic model. Such a model would 
have a softer equation of state, and would incorporate resonant decays. In order to change the ratio 7pb a i/7 p b a u one 
would have to spread the momentum balance around as a function of pt- Since 7 pbal has the more direct relation to 
the sum rule, it is less model-dependent. STAR certainly has the ability to determine 7 pba j, and has even analyzed 
the correlation differentially in terms of pt ■ 

Two main conclusions come out of the studies here. First, momentum conservation does seem to play a principle 
role in determining the combination of opposite- and same-sign correlations, 7 P bai and 7 p b a i- F° r a perfect acceptance, 
the rule seems unquestioned, except for the uncertainty in accounting for the fraction of momentum balance from 
unobserved particles (mainly neutral particles). Although the centrality dependence was not matched by the pion 
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cascade model, the size of the effect was very much of the right order, and for more central collisions the spread in 
A?7 was rather well reproduced. 

The second conclusion concerns the failure of the pion cascade to reproduce the strong centrality dependence 
of 7 P bai seen in the data. STAR's results suggest that the correlation is spread over a large range in rapidity for 
peripheral collisions. One intriguing explanation of the model's failure is that whereas the model evolution began at 
1 fm/c, for peripheral collisions it is possible that the elliptic flow was generated while the system was in a coherent 
state. If coherent structures, like color flux tubes or the color glass condensate, extend a unit of rapidity or more, 
momentum conservation might no longer be local (assumption #1 above) in r\. For this explanation to hold, much of 
the elliptic flow would have to be generated while the system were in such a coherent state. Whereas 1 fm/c represents 
a small fraction of the time over which flow develops in a central collision, the first fm/c should be more important 
for peripheral collisions, where the initial logarithmic density gradients are higher, and where the generation of flow 
finishes earlier. An even more speculative possibility is that in peripheral collisions elliptic flow might be generated at 
the initial point of contact between the nuclei. This might come from a preference for the individual nucleon-nucleon 
collisions to have their reaction planes coincide with the AA reaction plane, which could lead to initially anisotropic 
stress energy tensors with the momentum balance spread over multiple units of rapidity. However, our own quick 
analysis of Glauber models does not suggest such an effect exists. If early flow is indeed dominating the explanation, 
it would serve as motivation to reconsider some analyses of elliptic flow, based on hydrodynamic pictures with little 
initial flow. These analyses have often been used to determine the viscosity of the quark-gluon plasma [17 . 

Both STAR's experimental analysis and the model analysis presented here can be improved. The pion cascade 
model was adequate for illustrating the effects, but is too far from reality to make any serious quantitative conclusion 
from a comparison to data. The experimental analysis will be strengthened by better statistics, by a more detailed 
study of systematics (such as determining which method for v 2 is more appropriate for scaling 7 p bai) and by evaluating 
7pbai m addition to 7 P bai- It is difficult to determine how to fix many of the deficiencies of the model. More realistic 
models require hydrodynamic prescriptions for the intermediate stage as well as dynamic models of the field-dominated 
pre-equilibrium stages, thus incorporating the generation and propagation of momentum fluctuations in such models 
is necessary. Such advances are clearly challenging, but understanding the questions raised in the previous paragraphs 
would be valuable. This class of observables might ultimately provide insight into the state and nature of the pre- 
thcrmalized state, or might even assist in constraining bulk properties of quark-gluon plasma. 
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